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ABSTRACT 

We report on Chandra X-ray and Gemini-North near-infrared if-band observations 
of the eclipsing accretion-powered millisecond X-ray pulsar SWIFT J1749. 4-2807 in 
quiescence. Using the Chandra observation we derive a source position of Right Ascen- 
cion: 17:49:31.73 and Dcclination:-28:08:05.08. The position is accurate to 0'.'6 (90 per 
cent confidence). We find one source at a magnitude AT=18.44±0.03 with a position 
fully consistent with the accurate Chandra X-ray localisation and a second source at 
A"=19.2±0.1 that falls close to the edge of the error circle in the deep iC-band images. 
The presence of a few weaker sources as suggested by previous iJ-band observations 
presented in the literature cannot the ruled out. There is marginal evidence that the 
brighter of the these two sources is variable. Follow-up spectroscopy of this potential 
counterpart will show if this source is the true counterpart to SWIFT J1749. 4-2807. 
If so, baring the presence of complicating effects such as heating of the mass-donor 
star, it would allow for the mass of the neutron star to be measured through the 
measurement of periodic radial velocity variations. 
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stars: neutron 



1 INTRODUCTION 

Neutron stars provide a laboratory to test the behaviour of 
matter under physical conditions that are unattainable on 
Earth. The description of the relations between pressure and 
density of matter under the extreme conditions encountered 
in neutron stars (the equation of state; EoS) is one of the 
ultimate goals of the study of neutron stars already for 40 
years now. 

An excellent way to constrain proposed EoS is by find- 
ing massive neutron stars, as each of the EoS predicts a 
maximum mass for a neutron star above which it will col- 
lapse into a black hole (Lattimer & Prakash 20011. Neu- 
tron stars more massive than the canonical value of 1.4 
Mq can be found in X-ray binaries, as in those systems 
the neutron stars can accrete a substantial amount of mass 
from its binary companion. There is strong evidence that 
millisecond pulsars evolve out of low-mass X-ray binaries 



(LMXBs), see for instance Wijnands & van der Klis 
|Chakrabarty fc Morgan] ( |1998[ ) and [Archibald et al 



(1998). 



(2009) 



substantially is 1.5±0.2 Mq ( Ozel et al. 2012 see also Zhang 
et al.||2011[ ), there is a sample of millisecond pulsars where 



high neutron stars masses have been measured. For instance, 
the mass of the neutron star in FIRST J102347.6+ 003841 
has recently been determined to be 1.71±0.16 Mq (Deller 



et al.||2012[ ) and the millis econd pulsar PS R J0751+1807 
has a mass of 2.1±0.2 M ( |Nice et al.1|2005[ ). The most ac- 



curately measured mass of a massive neutron star is that in 
PSR J1614 -2230. It has a mass of 1.97+0.04 M Q ( jDemorest 
et al.|20T0 l. However, the maximum neutron star mass may 



be higher still. For instance, the best estimate of the Black 
Widow pulsar mass i s 2.40±0.12 M (PSR B1957+20 ;i _ 
Kerkwijk et al.|[2011 l, although the systematic uncertainty 



Whereas the average mass of neutron stars that accreted 



on this measurement could be larger than the statistical one. 

For dynamical mass measurements in LMXBs that host 
(millisecond) pulsars, pulse timing can directly constrain 
the orbit of the neutron star and thus its radial velocity 
semi-amplitude, K\. The radial velocity semi-amplitude of 
the mass donor star, K2, can be measured from optical or 
near-infrared observations. The ratio between K\ and K% 
yields the ratio between the mass of the companion star 
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and the neutron star, q. Alternatively, this q can be deter- 
mined from the rotational broadening of the stellar absorp- 
tion lines (v sin i) . This v sin i combined wit h Jfg gives q 



(1 + 9); 



0.49q 2 / 3 



(e.g. see 



jf^- - v- 1 - -r y-'o.69 2 /3 + i n{ i +9 i/3 ) v--5- |Horne et al. 
1986|). The system inclination can be determined through 



modelling of the multi-colour optical lightcurves or, in sys- 
tems with favorable viewing angles such as SWIFT J1749.4- 
2807 the X-ray eclipse duration can be used to accurately 
determine the inclination ( |Horne||1985[). As the inclination 
is constrained by the geometry ( Chanan et al.| | l97 6), mass 
measurements in eclipsing systems are independent of the 
modelling that lies behind inclinations derived from ellip- 
soidal variations. In most cases the system has to be in 
quiescence such that the accretion continuum is suppressed 
permitting the detection of photospheric absorption features 
from the mass donor. These absorption features allow one to 
measure K2 and vsini from optical or near-infrared spec- 
troscopic observations. The eclipse duration together with q 



gives an accurate measure of the inclination i ( Chanan et al 



1976 Home 1985j). q, i and K2 together give the neutron star 
mass. 



Recently, it was discovered that the transient 
SWIFT J1749. 4-2807 found in 2006 ( |Schady et aLl|2006| 
see also Wijnands et al.|2009 1 exhibits pulsations at 518 Hz 
( Altamirano et al. 2011 1. Furthermore, the X-ray light curve 
shows eclipses at the 8.8 hr orbital period ( jMarkwardt fc| 
Strohmaycr 20101. For the companion star to fill its Roche 



lobe in the 8.8 hours orbit its mean density would imply 
a spectral type of a G3 V-G5 V star if the companion is 
close to the lower main-sequence. Thus, this source holds 
the promise of allowing for a model-independent neutron 
star mass determination. From a type I X-ray burst the up- 
per limit to the distance to this source has been determined 
to be 6.7±1.3 kpc ( |Wijnands et al.|2009 l. 



The observed interstellar hydroge n column density Nj j 
towards the source is ~ 3 x 10 cm" " 2 ( |Ferrigno et al.|201l[ ). 
As it implies about 17 magnitudes of extinction in the V- 



band (Predehl & Schmitt 1995) it precludes the search of 



an optical counterpart. Given that the extinction in the K— 
band is only ~11% of that in the U-band ( |Cardelli et al.| 
1989 1 it allows for the detection of a near-infrared counter- 



part. D'Avanzo et al. (2011) searched for a near-infrared 



counterpart in the 1. 6 Swift error circle in the i/-band. 
Given the crowded field, the relatively large uncertainty 
on the Swift X-ray position and the deep observations, 
|D'Avanzo et al.| ( |2011[ ) found more than 40 potential coun- 
terparts. 



In this Paper we report on a Chandra observation of 
SWIFT J1749.4-2807 in quiescence. This observation pro- 
vides an accurate localisation. Using this accurate position 
we investigate our deep Gemini 7^-band observations ob- 
tained under excellent natural seeing conditions for the pres- 
ence of a counterpart to SWIFT J1749.4-2807. Throughout 
this paper we use the ephemeris of SWIFT J1749. 4-2807 



15" 



XMM 1 




.Chandra 
■ ■" ■ ■ 



Figure 1. Zoom-in around the Chandra ACIS-S3 position of 
SWIFT J1749. 4-2807 (indicated with the small red circle with 
a radius of 0"6). The blue circle indicates the Swift jD'Avanzo| 
|et al.|2011[ l position and the magenta circles represent the XMM- 
Newton pos itions of SWI FT J1749.4-2807 given by |Wijnands| 
|et al.||2009"l l (XMM 1) and |Degenaar et aljpull) (XMM 2). 



2 OBSERVATIONS, ANALYSIS AND RESULTS 
2.1 Chandra X-ray observation 

We observed SWIFT J1749. 4-2807 with the Chandra 
satellite using the back-illuminated S3 CCD-chip of 
the Advanced CCD Imaging Spectrometer (ACIS) detec- 
tor ( |Garmire 19971. The observation started on MJD 
56138.09427 (July 30, 2012, 02:15:44 UT) and the on-source 
exposure time was 24.4 ksec. The observation identification 
number is 13704. 

To mitigate potential pile-up on the off chance that the 
source was brighter than predicted we windowed the ACIS- 
S3 CCD such that only half the original CCD size is read out, 
providing a frame time of 1.54 s of which 0.04 s is used for the 
CCD read out. We have reprocessed and analysed the data 
using the CI AO 4-3 software employing the calibration files 
from the Calibration Database version 4.5. In our analysis 
we have selected events only if their energy falls in the 0.3-7 
keV range. All data have been used, as background flaring 
is very weak or absent. 

Since, by design, the source position falls near the op- 
tical axis of the telescope (see below), the size of the point 
spread function is smaller than the ACIS pixel size. There- 



fore, we follow the method of Li et al. ( 2004 1 implemented 



in the CIAO 4.3 tool ACIS_pr.OCESS_events to improve the 
image quality of the ACIS data. 

Using wavdetect we detect one X-ray source within 



the Swift (D'Avanzo et al. 2011 1 and XMM- Newton (Wi- 



given by Altamirano et al. (2011) 



jnands et al.||2009[ ) error circles of SWIFT J1749.4-2807 
(see Figure [lj. The source is detected 0.2' off-axis on the 
ACIS S3 CCD. The wavdetect J2000 source position in 
decimal degrees is: right ascencion (RA): 267.38220(2) and 
declination (Dec): -28.13474(1). The digit in-between brack- 
ets denotes the wavdetect 68 per cent confidence error on 
RA and Dec. The position in sexagesimal notation is: RA: 
17:49:31.73 Dec:-28:08:05.08. The errors given above are due 
to the error in the source localisation on the CCD alone. The 
absolute accuracy of the source position is thus determined 
by the uncertainty in the boresight of the Chandra satellite. 
For sources on the ACIS-S3 CCD this uncertainty is typi- 
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cally 0.6" at 90 per cent confidence. We tried to reduce this 
boresight uncertainty by investigating if one of the other de- 
tected X-ray sources can be safely associated with a source 
with a well known position but this was not the case. 

We have extracted source counts from a circular region 
of 2" radius centered on the source position. Similarly, we 
have used a circular region with a radius of 50" on a source- 
free region of the CCD to extract background counts. A 
point-source aperture correction was applied to the auxil- 
iary response file of the source. The net, background sub- 
tracted, detected source count rate is 2.9xl0~ 3 counts s _1 . 
In total after correcting for the number of expected back- 
ground photons, 61 source photons have been detected in 
the energy range between 0.3-7 keV. The number of back- 
ground photons in the same area in 0.3-7 keV was less than 
1. 

Using xspec version 12.6.0q ( |Arnaud|1996"l ) we have fit- 
ted the spectra of SWIFT J1749.4-2807 using Cash statis- 
tics |Cash 1979) modified to account for the subtraction of 
background counts, the so called W-statistic^] 

Given that we only detected a low number of counts 
we decided to fix the Nh to 3xl0 22 cm -2 in our spectral 
fits, consistent with the values determined by |Ferrigno et alT| 
(20111 and Degenaar et al. ( 2012 1 . We have used an ab- 



sorbed powerlaw model fit function to describe the data (the 
model pegpwrlw in xspec). We visually inspected a plot of 
the count rate as a function of heliocentric time to search for 
the presence of the eclipse. The Chandra observed count rate 
is consistent with zero at the predicted heliocentric time of 
the eclipse. The eclipse duration is consistent with the dura- 
tion determined by |Markwardt fc Strohmayer] ( |2010| ) when 
the source was in outburst. To account for the fact that the 



source is eclipsed during 2 ksec (Markwardt & Strohmayer 
|2010[ ), which is about 8 per cent of the total exposure time, 
we decrease the effective exposure by 2 ksec when calculating 
the source flux. We find a best-fit power law index of 0.6±0.4 
which gives an unabsorbed X-ray flux of (0.8 ±0.1) x 10~ 13 
erg cm -2 s -1 in the 0.3-7 keV band and 1.4 x 10~ 13 erg 
cm -2 s _1 in the 0.5-10 keV band. 



2.2 Gemini near— infrared observations 

if -band imaging of the field containing SWIFT J1749.4- 
2807 was performed with the Gemini-North telescope on the 
nights of June 15 and 20, 2012 under observing program GN- 
2012A-Q-67. The Near Infrared Imaging and Spectrometer 
(NIRI; jHodapp et aT | [2003] ) instrument was used with the 



f/14 camera to yield a field of view of 51" x 51" with a plate 
scale of 0'.'05 pixel" 1 . 

The field of SWIFT J1749.4-2807 is crowded and care 
was taken to avoid as much as possible image artifacts due 
to saturated stars. In particular, we exclude from the NIRI 
FOV the K ~ 8.6 star 2MASS J17493218-2807563. There- 
fore, we both rotated the instrument 45 degrees and used a 
(p,g=x,0)-dithering pattern. Integration times were 10 s in 
other to mitigate saturation effects from other bright stars 
in the field of view. Thus the observing sequences consisted 
of a 5-point dither pattern with offset steps p of 0, 10, -5, 5 
and 15" and we added three 10 s exposures at each dither 




Figure 2. Gemmi-NIRI K band finder chart obtained from com- 
bining those images with seeing better than 0"31 on June 15, 
2012. The 0"6 Chandra 90 per cent confidence error circle on 
the X-ray source position is indicated (circle). We mark source 1 
and 2, where the position of source 1 is consistent with the X-ray 
position. Source 1 shows marginal evidence for variability (at the 
4 <t level) between the observations of June 15 and June 20, 2012. 



position. This pattern was repeated twenty one times yield- 
ing a total of 52.5 min on source. The image quality varied 
during the observations due to the weather conditions. We 
measure a point-spread function full-width at half maximum 
(FWHM) between 0'.'22-0'.'36 and 0"23-0"58 with a median 
value of 0'.'29±0'.'03 and 0'.'40±0'.'07 on June 15 and 20, re- 
spectively. 

The data were processed using both python scripts and 
IRAlP] tasks developed and provided by the Gemini Obser- 
vatory. Note that sky frames were not subtracted from the 
data. The observed field is too crowded to allow us to build 
a useful sky frame from the dithered images and there were 
no nearby regions on the sky suitable to measure the local 
sky background. 

For the analysis, we generated images by stacking all 
the reduced frames using imcoadd. In this way we obtained 
for the night of June 15 seven combined images of 150 s 
on source each. We also combined the images with seeing 
<0"31 of the first night into one average image. The night 
of June 20 was of less good image quality (see above) and a 
single combined frame was build totalling 1770 s on source. 

We improved the astrometry of the best-seeing frame 
using the position of six 2MASS stars that were well sepa- 
rated from other stars and which did not saturate the Gem- 
mi NIRI detector. A fit of pixel scale in x and y, and ori- 
entation provide an root-mean-square (rms) uncertainty of 
0"03 for a pixel scale of 0"05 per pixel in x and y and an 
orientation of 44.9°. 

In Fig. [2] we show the finding chart of the field of 
SWIFT J1749.4-2807 using the best seeing image of June 
15. Overplotted is the 90 per cent confidence error circle on 
the Chandra X-ray position of the source. A clear source 
is detected near the centre of the error circle (star 1) and 
towards the south-east corner just outside the formal 90 
per cent confidence error circle, lies another fainter poten- 
tial candidate counterpart (star 2). 



1 see http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/ 



2 I RAF is distributed by the National Optical Astronomy Obser- 
vatories 
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We performed aperture photometry on the images us- 
ing daophot in iraf to compute the instrumental magni- 
tudes for the objects of interest. Absolute flux calibration of 
the images was obtained selecting isolated stars with magni- 
tudes reported in the UKIDSS-DR6 Galactic Plane Survey 
(UGPS; |Lucaa et al.|2008) . Both the NIRI and the UKIDSS 
photometric systems are close representations of the Mauna 
Kea Observatory (MKO) standard system. Differential pho- 
tometry was obtained on each of the averaged images to de- 
rive the flux variability of the sources as a function of time. 
The photometric results given here are with respect to the 
star UGPS J176971. 32-288889. 6 (K = 14. 551 ± 0.022). The 
scatter in the light curves build for this UGPS star using 
other bright and isolated comparison stars shows that the 
differential photometry is accurate to ^ 1 per cent. 

Photometry of the two candidate counterparts on the 
combined best seeing image of June 15, 2012, yields a mag- 
nitude of Ki = 18.44 ±0.03 and K 2 = 19.20 ±0.06 for star 1 
and star 2, respectively. The magnitude during the observa- 
tion on June 20, 2012 is Ki =18.75±0.07 and K 2 = 19.3±0.1 
mag. So, with an average magnitude of K\ = 18.44 ± 0.03 
on June 15, and K x = 18.75±0.07 on June 20, 2012, there is 
evidence at the 4 a level that this source is variable. 



3 DISCUSSION 

Using Chandra we obtained a 24.4 ksec exposure of the 
eclipsing msec, pulsar SWIFT J1749. 4-2807 in quiescence. 
We detected the source at an unabsorbed flux of 2 x 10 -13 
erg -1 cm -2 s -1 in the 0.5-10 keV band and we provide a 
position accurate to 0'.'6 at 90 per cent confidence. At a dis- 



tance upper limit of 6.7±1.3 kpc (Wijnands et al 



2009[ ), 

this flux corresponds to an X-ray luminosity of lxlO'" erg 
s _1 in the 0.5-10 keV band. Owing to the uncertainty in 
the high hydrogen column density towards the source this 
number is uncertain, but the value is consistent with that 
found before from XMM-TVeuiton observations in quiescence 
(Wijnands et al.||2009] |Degenaar et al.||2012[). The derived 



quiescent X-ray luminosity is also in agreement with the 
rough correlation between quiescent X-ray luminosity and 
the orbital period, although the source is on the bright end 
of the observed correlation (a plot highlighting this corre- 
lation can be found in |Rea et aij]2011| >. The best-fit power 
law index of the X-ray spectrum is 0.6 ± 0.4, which is steep 
for quiescent neutron stars, however, the uncertainty on this 
number is currently such that it could be made consistent 
with more canonical values of 1.5-2 particularly if Nh is 
somewhat higher than the 3xl0 22 cm -2 that was used in 



the fit (see for instance the powerlaw fit in Degenaar et al 



2012). 



In addition to the Chandra observation, we obtained 
Gemini _K-band images under excellent natural seeing con- 
ditions of the field of SWIFT J1749.4-2807. The position 
of one point source is consistent with the X-ray location 
and a second is bordering the Chandra error circle. There is 
marginal evidence for the brighter of the two to be variable. 
The observed interstellar hydrogen column density, Nh, to- 
wards SWIFT J1749. 4-2807 



3 x 



et al. 



mag 



I Ferrigno 



2011 1. This yields an extinction in the K-b&nd of 1.9 
Cardelli et al.|1989| ). However, as A V /E(B - V) = R v 
is known to be closer to 2.5 than 3.1 for fields towards the 



bulge ( Nataf et al.|2012| ) this implies that Ak would be 1.7 
mag for iiy=2.5. 

For the companion star to fill its Roche lobe in the 8.8 
hours orbit its mean density would imply a spectral type of 
a G3 V-G5 V star if the companion is close to the lower 
main-sequence. Note that the spectral type of G5 V that 
we derive here is consistent with the broad mass range for 
the companion star (0.45-0.8 Mq) that is derived by Mark 



wardt & Strohmayer (20101 and Altamirano et al. (20111 



as the mass donors in LMXBs and CVs are often under- 
massive for their spectral type (iKnigge et al. 20111. For a 



G5 V donor star, the absolute magnitude in the V- and in 
the if-band is My =5.1 and Mk = 3.5, respectively (see 
Cox||2000 and Tbkunaga 2000, respectively). At a distance 
upper limit of 6.7 kpc the distance modulus, DM, is 14.1. 
Allowing for 1.7 magnitudes of extinction in the if-band 
will mean that the companion star has an apparent magni- 
tude of mx=DM±1.7±M K yielding m K = 19.3 for the G5 V 
star. The observed K band magnitude of ~ 18.6 is consis- 
tent with this if the source is more nearby than the 6.7 kpc 
upper limit (i.e. the source has to be at ~5 kpc) and/or if 
the extinction in the A'-band is lower than 1.7 magnitudes. 
The latter can be the case if part of the observed Nh is due 
to ionised material for instance material local to the source. 

Previous to our work, |D'Avanzo et al.| ( |2011[ ) obtained 
//-band adaptive optics observations with the Very Large 
Telescope. Of the 41 sources that those authors detected 
inside the error circle of the Swift X-ray position, five fall 
inside the Chandra error circle; six if we include our source 
2. Two of these six sources coincide with the position of the 
two potential counterparts that we found. In Table [l]we list 
the properties found by D'Avanzo et al. (2011 1 for these six 
sources. 

Our star 1 at RA: 17:49:31.73 Dec: -28:08:05.09 would 
correspond to source number 8 in table 2 of |D'Avanzo et aT7| 
( |2011[ ) (see Table[l]). Our star 2 which is at RA: 17:49:31.768 
Dec: -28:08:05.41 would correspond to their entry number 
24. The other four sources inside the Chandra error circle 
are not significantly detected in our Gemini data, although 
there is some evidence for additional flux to the south of 
star 1, which could be caused by either source 25 or source 
27, or by both sources together. 

The extinction towards SWIFT J1749.4-2807 in the 
//-band is ~1.2 mag more than that in the /f-band. The 
(H — K)o colour of A-M5 stars is between and 0.3 mag 
( jTokun aga 2000j. Therefore, the magnitudes of the two 
sources that we find should be, at most, 1.5 magnitude 
brighter than the //-band magnitudes reported in |D'Avanzo| 
|et al.| ( |2011[ ) for these two sources. However, comparing our 
observed K band magnitudes with the //-band magnitudes 
reported by D'Avanzo et al. (2011 \ (see Table[l]) we see that 
the observed H — K colour is approximately 2.4 magnitudes 
for star 1 and 2.8 magnitudes for star 2. We investigated 
the potential reason for this peculiar colour by comparing 
our magnitudes as well as the //-band magnitudes derived 



by D'Avanzo et al. (20111 with the magnitudes of 2MASS 



and UKIDSS stars in the images. The //-band magnitude 
of the star at RA : 17:49 :31.88 Dec: -28:08:03.3 reported by 
Avanzo et al. (20111 (entry number 1 in their paper) 
is 16.26, however it has //=15.01±0.05 according to the 
UKIDDS database, this is nearly 1.3 magnitudes brighter. 
In addition, we looked up the H and K magnitudes reported 
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Table 1. Position and //-band magnitudes from [jTAvanzo et al.| j2011| for the sources close to the Chandra 90 per cent confidence 
error circle. The offset of the position of each of the stars with respect to the centre of the Chandra position as well as the /sf-band 
magnitudes from this work are given. Note that the //-band magnitudes from |D'Avanzo et al.||2011| are off; approximately —1.3 has to 
be added to the reported magnitudes (see Discussion). 



RA 


Dec 


Src # 


//-band mag 


//-band mag 


Offset 


K— band 








D'Avanzo et al. 2011 


Aug 30, 2010 


Aug 31, 2011 


a 


mag 




17:49:31.730 


-28:08:04.50 


32 


22.75±0.11 


22.66±0.09 


0.58 






17:49:31.76 


-28:08:05.39 


41 


23.99±0.35 


22.68±0.09 


0.53 






17:49:31.720 


-28:08:05.50 


25 


22.22±0.07 


21.86±0.07 


0.44 






17:49:31.730 


-28:08:05.40 


27 


22.54 ±0.08 


22.38±0.08 


0.32 






17:49:31.72 


-28:08:05.1 


8 


20.81±0.05 


20.74±0.05 


0.11 


18.44±0.03 (src 


1) 


17:49:31.77 


-28:08:05.4 


24 


22.20±0.08 


21.94±0.12 


0.64 


19.20±0.06 (src 


2) 



stars are the CO band head absorption bands (e.g. Greiner 




et al. 2001b Greiner et al. 2001a I. If confirmed, source 1 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Orbital phase (1.5 cycles shown) 

Figure 3. The expected ellipsoidal variations for a Roche lobe 
filling G5V star with in decreasing order of amplitude a mass ratio 
of q = 0.3, 0.5, 0.7. The magnitudes for source 1 on the night of 
June 15, 2012 (near phase 0.77) and June 20, 2012 (near phase 
0.36) are overplotted, showing that the observed variability has 
the right sign in the sense that the source is brighter near phase 
0.77 than near phase 0.36. Changes in the predicted ellipsoidal 
modulations for a K0 V star are marginal with respect to those 
of the G5 V star shown. 



in the VISTA (Visible and Infrared Survey Telescope for As- 
tronomy) variable sources in the Via Lactea (VVV) survey 
( Minniti et al.|2010|[Catelan et al.|2011[ ) and the bright star 
has £Z=14.9, so close to the value provided by UKIDSS. 
We conclude that the photometric calibration presented in 
D'Avanzo et al. (20111 is off by 1.3 magnitudes. Indeed, 



correcting for this photometric offset would imply observed 
H — K colours of our stars 1 and 2 in line with what is 
expected given the extinction towards the source. The alter- 
native scenarios that could potentially explain the observed 
H — K colours, such as a significantly higher extinction or 
a more shallow extinction law (a large value of Rv) are not 
likely nor necessary to explain the H — K colours. 

As a next step, spectroscopic follow-up is suggested to 
confirm/reject source 1 as the counterpart. Spectra should 
show typical late type star absorption features that should 
move as a function of the orbital phase of SWIFT J1749.4- 
2807. Typical features in the if -band spectra of late type 



is within reach of current instrumentation for time-resolved 
dynamical studies. If source 1 is found not to be the coun- 
terpart, potentially, one can search for the right counter- 
part by using adaptive optics observations under photomet- 
ric conditions. However, given that the maximum amplitude 
(peak-to-peak) of ellipsoidal variations is 0.2 magnitudes 
for a mass ratio of q =0.3 this is challenging. In Figure[3]we 
plot the magnitudes of source 1 over the expected ellipsoidal 
variations for a Roche lobe filling companion star of spec- 
tral type G5V for three different values of q. With decreas- 
ing amplitude of the modulations we plot q — 0.3, 0.5, 0.7. 
The figure shows that the observed variability has the right 
sign given that the source is brighter near phase 0.77 than 
near phase 0.36. For smaller mass ratios the peak-to-peak 
amplitude increases, which given that the two data points 
lie above and below the curves, could hint at a more ex- 
treme mass ratio. Such a mass ratio could be accommodated 
by a neutron star more massive than 1.4 Mq and a mass 
donor star under-massive for its spectral type. Furthermore, 
if light from a residual accretion disc is present it will reduce 
the apparent amplitude of the ellipsoidal modulations. Thus, 
adding a significant amount of light of an accretion disc to 
the model can only be made in agreement with the two data 
points for a (very) small mass ratio q. Potentially, if the neu- 
tron star heats the companion star one can find variations 
due to changes in the aspect of the companion star (see 
for instanc e |Homer et al.||2001| |Burderi et aTj]2003[ |Jonker| 
|et al.|2008| |D'Avanzo et al.|2009 |. However, such variations 
would possibly be larger than the ellipsoidal variations and 
they would produce a bright if -band magnitude at orbital 
phase of 0.5 and this is not seen in the current data (see 
Figure|3|. Given that the orbital period of SWIFT J1749.4- 
2807 is larger than that of the accretion powered millisecond 
pulsars where this effect has been found before, the ampli- 
tude of the effect will be reduced. Finally, flickering often 
observed in optical and near-infrared light of X-ray binaries 
in quiescence due to for instance short-term variations of 
the accretion disc can also explain the observed variability 
(e-g- 
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